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Abstract We present an algorithm for computing the set of all coset leaders of a 
binary code 'lo C F". The method is adapted from some of the techniques related 
to the computation of Grobner representations associated with codes. The algorithm 
provides a Grobner representation of the binary code and the set of coset leaders 
CL('^). Its efficiency stands of the fact that its complexity is linear on the number of 
elements of CL('^), which is smaller than exhaustive search in Fj. 

Keywords Binary codes • Cosets leaders ■ Grobner representations. 
1 Introduction 

The error-correction problem in coding theory addresses given a received word recov- 
ering the codeword closest to it with respect to the Hamming distance. This previous 
statement is the usual formulation of the Complete Decoding Problem (CDP). The 
f-bounded distance decoding (f-BDD) algorithms determine a codeword (if such a 
word exists) which is at distance less or equal to f to the received word. If t is the 
covering radius of the code then the bounded distance decoding problem is the same 
as CDP. In the CDP of a linear code of length « , ^ C FJ^ those errors that can be 
corrected are just the coset leaders, which are vectors of smallest weight in the cosets 
F^/*^. When there is more than one leader in a coset there is more than one choice 
for the error. Therefore the following problem naturally follows known as the coset 
weights proble}n{CWP): 
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Input: A binary r xn matrix H, a vector s S and a non-negative integer t. 
Problem: Does a binary vector e 6 of Hamming weight at most t exist such that 

Recall that H can be seen as the parity check matrix of a binary code and s the syn- 
drome of a received word, thus the knowledge of e would solve the f-BDD problem. 
Unfortunatelly the hope of finding an efficient f-BDD algorithm is very bleak since it 
was proven that the CWP is NP-complete 01 ■ Thus the computation of all the coset 
leaders is also NP-complete. The study of the set of coset leaders is also related to 
the study of the set of minimal codewords, which have been used in the Maximum 
Likelihood Decoding Analysis |2J and which are also related to the minimal access 
structure of secret sharing schemes ifTSl . Furthermore, the computation of all coset 
leaders of a code allows to know more about its internal structure Q- 

All problems mentioned before are considered to be hard computational problems 
(see for example EHl) even if preprocesing is allowed ifTOl . However, taking into 
account the nature of the problem, to develop an algorithm for computing the set of 
all coset leaders of a binary code, in the vector space of 2" vectors, by generating 
a number of vectors close to the cardinality of this set may be quite efficient. This 
is our purpose extending some results on Grobner representations for binary codes. 
For previous results and applications of Grobner representations for linear codes see 
SEMI, for a summary of the whole material we refer the reader to IS]. We extend 
some settings of previous work in order to obtain the set of all coset leaders keeping 
record of the additive structure in the cosets. The presentation of the paper is done in 
a "Grobner bases"-free context. 

The outline of the paper is as follows. In Section|2]we review some of the stan- 
dard facts on binary linear codes and their Grobner representation. Section [3] gives 
a concise presentation of the results in this paper. Definition [3] corresponds to the 
construction of List, the main object that is used in the algorithm proposed, while 
Theorem |2] guarantees that all coset leaders will belong to List. In this section it is 
presented the algorithm CLBC for computing the set of all coset leaders. At the end 
of the section we show an example to illustrate a computational approach applied to a 
binary linear code with 64 cosets and 118 coset leaders. In Section|4]we discuss some 
complexity issues. Finally some conclusion are given which include further research. 



2 Preliminaries 

Let F2 be the finite field with 2 elements and Fj be the F2-vector space of dimension 
n. We will call the vectors in Fj words. A linear code '^^ of dimension k and length 
n is the image of an injective linear mapping L : Fj Fj, where k <n, i.e. = 
L{¥\). From now on, we will use the term code to mean binary linear code. The 
elements in '^i' are called codewords. For a word y G F" the support of y is supp(y) = 
{i G {1, . . . ,«} I y; ^ 0} and the Hamming weight of y is given by weight(y) the 
cardinal of supp(y). The Hamming distance between two words Ci,C2 is d(ci,C2) = 
weight(ci — C2) and the minimum distance li of a code is the minimum weight among 
all the non-zero codewords. It is well known that CDP has a unique solution for those 
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vectors in ,t) = {y £ \ 3c<^'^ s.t. c/(c,y) < [^] } where [•] is the greatest 
integer function. 

Definition 1 The words of minimum Hamming weight in the cosets of F"/*^ are 
called coset leaders. 

Cosets corresponding to B{'^^,t) have a unique leader however in general outside 
B{'^,t) there may be also cosets with a unique leader, i.e. those cosets with only one 
leader could be more than |B('^,f)|. Let CL('^) denote the set of coset leaders of the 
code and CL(y) the subset of coset leaders corresponding to the coset +y. Let e, 
; = 1, . . . ,« be the i-th vector of the canonical basis X = {ei, . . . , e,,} C Fj and the 
zero vector of Wj. The following theorem gives some relations between cosets lfT2l 
Corollary 11.7.7]. 

Tlieorem 1 Let w G CL('^^), such that w = W] +ei for some wi S F" and ; g { 1 , . . . , «}. 
Then, wi G CL(wi). 

Definition 2 A Grobner representation of F"/*^ OH] is a pair N,0 where N is a 
transversal of Fj/*^ such that G and for each n G N \ {0} there exists e,-, / G 
such that n = n' + e,- and n' G N, and (/) : N x {e,}JL[ N be the function 
that maps each pair (n,e, ) to the element of N which belongs to the coset of n + e,. 

3 Computing tlie set of coset feaders 

A key point of the algorithms for computing Grobner representations is the construc- 
tion of an object we called List which is an ordered set of elements of Fj w.rt. a 
linear order ^ defined as follows: w ^ v if weight(w) < weight(v) or weight(w) = 
weight(v) and weight(w) weight(v), where -<i is any admissible order on Fj in 
the sense given in ||3] p. 167]. We will call such kind of orders as weight compatible 
ordering s. 

Definition 3 (Construction of List) . Let List be the ordered structure given by the 
following axioms 

L OGList. 

2. If V G List, letN(v) =min^{w | w G Listn + v)}. 

3. If V G List is such that weight(v) = weight(N(v)), then {v + e; : ; ^ supp(v), / G 

CList. 

Remark 1 The set N is the subset of List such that N(v) is the least element in Listfl 
(^ + v), i.e. N(v) G ''^ + v andN(v) =min^{w | w G Listn C'^ + v)}. 

Remark 2 We start List with and we will see that, whenever condition 3. holds for 
V, the vector v will be a coset leader of + \. Then, we insert v + e, to List, for 
i — 1 , . . . , « and / ^ supp(v) (see Theorem[T]i. It is not necessary to introduce v + e,, 
for / G supp(v), because in this case v + e, -< v and then v + e, has been already 
considered in List since it is an ordered structure. 
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Next theorem states that List in Definition|3]includes the set of coset leaders. 
Theorem 2 Lef w G CL{'^^), then w £ List. 

Proof We will proceed by Noetherian Induction on the words of with the ordering 
since -< has the property of being a well- founded ordering or Noetherian ordering 
(i.e., any descending chain of words is finite) pi Theorem 4.62, p. 168]. The word 
G CL('i') and by definition it belongs to List, let w G CL('^) \ {0}. Assume the 
property valid for any word less than w with respect to -<, i.e. u G List provided 
u G CL('^^) and u ^ w. Taking ; G supp(w), we write w = u + e,- where u G Fj, 
then by Theorem [T] u G CL('^^). In addition, supp(u) — supp(w) — 1, thus u w. 
Therefore, by applying the induction principle we have u G List. If u is a coset leader 
belonging to List it is clear that weight(u) = weight(N(u)). As a consequence, by 3. 
in Definition[3] w = u + e, G List. □ 



Algorithm 1 (CLBC) 

Input: A parity check matrix of a binary code 

Output: CL('^), 0: The set of all coset leaders and the function Matphi. 

1: List ^ [0], N ^ 0, r ^ 0, CL{'^) ^ [] 

2: while List 7^ do 

3: T <— NextTerm[List], s <^ xH 

4: y Member[s, {si, . . . ,s, }] 

5: it j false then 

6: for k such that T = t' + e^; with t' G N do (t' , e^; ) 

7: if weight(T) — weight(Ty) then 

8: CL(^)[Tj] ^ CL{'^)[tj] U {t} 

9: List ^ InsertNext[T, List] 

10: end if 

11: else 

12: r ^ r+ 1, Sr ^ S, T;- ^ T, N ^NU{t,.} 

13: CL{'^)[Tr] ^ CLC^) [t,] U {t} 

14: List <— InsertNext[Tr,List] 

15: for k such that T,. = t' + e^; with t' G N do 

16: (j){T',ek) 

17: 0(T;-,eA.) ^ t' 

18: end for 
19: end if 
20: end while 

21: return CL^^), (j) 

Where 

1. InsertNext[T,List] Inserts all the sums T+e^ in List, where k ^ supp(T), and keeps 
List in increasing order w.r.t. the ordering 

2. NextTerm[List] returns the first element from List and deletes it from this set. 

3. Member [oZ? j, G] returns the position j of ob j in G if obj E G and false otherwise. 
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Theorem 3 CLBC computes the set of coset leaders of a given binary code and its 
corresponding Matphi. 

Proof Note that when an element T is taken out from List by NextTerm the elements 
to be inserted in List by InsertNext are of the form T + e^;, where k ^ supp(T). As 
a consequence, T + e^ >- T. Then all elements generated by CLBC in List, after T is 
taken out, are greater than T. Therefore, when T is the first element in List in Step 3, 
all elements of List that shall be analyzed by CLBC are greater than T. 

Let us prove that the procedure generates List according to Definition[3] First, by 
Step 1, e List. Let T = NextTerm [List] in Step 3, in this step the syndrome s of T is 
computed. Thus we have two cases regarding the result of Step 4, namely 

1. Assume j = "false", thus Steps 5 and 11 guaranty us to find N(t) as the least 
element in List having the same syndrome as T. When we are in this case (t = 
N(t)), in Step 14 it is performed (ii) of Definition [3] 

2. On the other hand, assume j ^ "false" (an element N(t) = T, have been already 
computed), if we have weight(T) = weight(T,), in Step 9 it is performed (ii) of 
Definition |3] 

Therefore CLBC constructs the object List following Definition[3] By Theorem|2] the 
set of coset leaders is a subset of List. Then Steps 8 and 13 assure the computation 
of this set. The procedure computes the Matphi structure as a direct consequence of 
Steps 6, 16, 17. The reason for including Step 17 is that in this case 

Tr + ei = t' -< T,- 

and those elements Xr + ^k, where k G supp(Tr), are not inserted in List since t' has 
been already considered when T, is computed as a new element of N. In addition, 
since -< is admissible, all subwords of are also least elements of their cosets ac- 
cording to ^, so t' = N(t') (note t' is a coset leader and by Theorem|2]it belongs to 
List). 

We have proved that CLBC guarantees the outputs we are claiming. Termination 
is a consequence of the fact that the cardinal of the set of elements belonging to List 
is least than n|CL('^#')|. Then, to a certain extent (when the set of coset leaders has 
been computed) no more elements are inserted in List in Steps 9 and 14. Therefore, 
the list get empty, and by Step 2, CLBC terminates. □ 

Remark 3 

1. We introduce in this paper steps 7, 8 and 9. In our predecessor algorithms where 
it was not necessary to introduce them ||6l|7][8l. These steps, as well as a modified 
definition of list, are forced by the fact that all coset leaders are incorporated. 

2. In Steps from 14 to 17, the pairs {Tr,^k), where k G supp(Tr), are not necessary 
to compute the coset leaders but for computing the structure Matphi and could be 
removed if we are only interested in the coset leaders. 

Example 1 Consider the [10,4, 1] code with parity check matrix H 
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1110000100 ■ 
1111000010 

yiiiioooooi/ 

We use GAP 4.12 fn\ and the GAP's package GUAVA 3.10 for Coding The- 
ory. We have built in this framework a collection of programs we call GBLA_LC 
"Grobner Bases by Linear Algebra and Linear Codes" Q- In particular, we have run 
the function "CLBC" of GBLA_LC (Coset Leaders of Binary Codes), it gives a list 
of three objects as an output, the first one is the set of coset leaders, the second one 
the function Matphi, and the third one the error correcting capability of the code. The 
complete set of coset leaders is the list below with 64 components and each compo- 
nent corresponds to a coset with its cosets leaders, the set N is composed by the first 
elements of each component. We have indicated with arrows some places in the list 
below, which we are going to use during the example. The elements in the list CL('if ) 
are 

[ [1] , [ei] , [e2] , [es] , [sa] , [es] , [eg] , [ey] , [eg] , [eg] , [eio] , 

[ei +62,^ es + 66 H Jei +63,65 + 67], [ei +64,65+68], 

[61 +65,62 +66,63 +67,64 + 68], [61 +66,62 +65], [61 +67,63+65], 
[6 1+68, ^64 + 65^], [6 1+69], [6 1+610], [62 + 63, 66 + 67], 

[62 + 64, 66 + 68] , [62 + 67, 63 + 66] , [62 + 68, ^ 64 + 66 ^] , [62 + 69] , [62 + 610] , 

[63 + 64, 67 + 68] , [63 + 68 , 64 + 67] , [63 + 69] , [63 + eio] , 

[64 + 69] , [64 + 610] , [65 + 69] , [65 + 610] , [66 + 69] , [66 + 610] , 

[67 + 69] , [67 + 6 10] , [68 + 69] , [68 + 6 1 0] , [69 + 6 1 0] , 

[61 +62 + 63,61+66 + 67,62+65+67,63+65+66], 

[61 +62 + 64,61 +66 + 68,62 + 65+68,64+65+66] 
[61 +62 + 67,61+63+66,62+63+65,65+66+67], 
[61 +62 + 68,61 +64 + 66,62+64 + 65,65+66+68], 
[61 +62+69, 65+66 + 69], [61 +62 + 610,65+66 + 610], 
[61 +63+64,61 +67+68,63+65+68,64+65+67], 
[61 +63+68,61+64 + 67,63+64 + 65,65+67+68], 

[61 +63 +69,65 +67 +69], [61 +63+610,65+67+610], [61 +64+69,65 +68+69], 
[61 + 64 + 610, 65 + 68 + 610] , [61 + 65 + 69, 62 + 66 + 69, 63 + 67 + 69, 64 + 68 + £9] , 
[61 + 65 + 610, 62 + 66 + 610, 63 + 67 + eiO, 64 + 68 + 610] , [61 + 66 + 69, 62 + 65 + 69] , 
[61 +66 + 610,62+65 +610], [61 +67 +69,63 +65 +69], [61 +67+610,63+65+610], 
[61 + 68 + 69, 64 + 65 + 69] , [61 + 68 + 610, 64 + 65 + 610] , [61 + 69 + 610] , 
[62 + 63 + 68 , 62 + 64 + 67 , 63 + 64 + 66, 66 + 67 + 68] , [65 + 69 + 610] ] 

Note that y = 64 + 65 + 66 in CL{^^)n(, = [ei + 62 + 64,61 + 66 + 68,62 + 65 + 
68 , 64 + 65 + 66] (pointed by arrows) is a coset leader such that none subword of y 
is in N (see the previous elements pointed by arrows). This shows the importance 
of considering all coset leaders in 3. of Definition [3] and not only the coset leaders 
belonging to N like in previous works. 
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The algorithm could be adapted without incrementing the complexity to get more 
information like the Covering radius and the Newton radius. In this case by analyzing 
the last element CL('^)64 = [es +eg +eio] we have a coset of highest weight which 
also contains only one leader; therefore, the Covering radius and the Newton radius 
are equal to 3. Moreover, for computing these parameters the algorithm do not need 
to run until the very end. 

As we state before, the set N would be enough to compute the Weight Distribution 
of the Coset Leaders WDCL = (oCq, . . . , OC,,) where a,- is the number of cosets with 
coset leaders of weight /, ; = 1 of the code. We can provide also an object which 
gives more information about the structure of the code, that would be the numbers of 
coset leaders in each coset (#(CL)). 

WDCL = [1,10,30,23,0,0,0,0,0,0,0]. 

#(CL) = [1,1, 1,1, 1,1, 1,1, 1,1, 2,1, 1,1, 2,2,4,1, 1,1, 1,2,2,2,2,2,2,4,2,2,2,4,4, 
1,1,1,4,4,4,2,2,2,2,2,4,4,1,1,1,2,2,2,2,4,1,1,1,2,2,2,1,1,1,1]. 

It is also interesting to note that there are more cosets with one leader (19) out the 
cosets corresponding to B(^^,r) than cosets corresponding to B(^,f) 

|{CL(y)|yeBr^,f)}| = ll. 

Therefore, there are 30 of the 64 cosets where the CDP has a unique solution. 



4 Complexity Analysis 

For a detailed complexity analysis and some useful considerations from the compu- 
tational point of view we refer the reader to f7\. Section 6 of that paper is devoted 
to discuss in details about the computational complexity and space complexity of 
the setting for computing Grobner basis representation for binary codes and general 
linear codes. The method is also compared with other existent methods for similar 
purposes. In the case of this paper, the difference is that we work out the set of all 
coset leaders and not only a set of canonical forms. Next theorem shows an upper 
bound for the number of iterations that will perform CLBC. 

Theorem 4 CLBC computes the set of coset leaders of a given binary code of 
length n after at most ii|CL('^) | iterations. 

Proof Note that the number of iterations is exactly the size of List. In the proof of 
Theorem|3]was shown that the algorithm follows this definition to construct the object 
List. It is clear that the size of List is bounded by n|CL('^)|, note that we can write 
List, as a set as follows 

List = {w + e,- I w e CL(^) and G {!,...,«}}. 

□ 
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Remark4 1. By the proof above we require a memory space of (9(«|CL(^)|). We 
assume that for computing the set of coset leaders it is required at least 0{ |CL('^^) | ) ; 
therefore, CLBC is near the optimal case of memory requirements. 

2. CLBC generates at most «|CL('^)| words from Fj to compute the set of coset 
leaders. An algorithm for computing this set needs to generate from F" at least a 
subset formed by all coset leaders, i.e. |CL('^^)|; therefore, the algorithm is near 
the optimal case of computational complexity. 

5 Conclusion 

The Algorithm CLBC is formulated in this paper, which turns out to be quite effi- 
cient for computing all coset leaders of a binary code from memory requirements and 
computational complexity view. Although, as it is expected, the complexity of the 
algorithm is exponential (in the number of check positions). 

The difference of CLBC with its predecessors rely also in the computation of 
all coset leaders instead of a set of representative leaders for the cosets, however, it 
supports the computation of the function Matphi. We remark that the computation 
of Matphi is not necessary at all for the main goal of computing the coset leaders. 
We have kept this resource in the algorithm because this structure provides some 
computational advantages (see 01381) and, although the algorithm will be clearly 
faster without computing Matphi, the nature of the computational complexity and 
space complexity will remain the same. 

Unfortunately, a generalization to non-binary linear codes is not trivial from this 
work. The main reason seems to be that the solution based on Hamming weight com- 
patible orderings will not continue being possible; the error vector ordering we have 
used in the general approach Q is not a total ordering, although it allowed to set 
up the computational environment in order to compute Grobner representations for 
Unear codes. 
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